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The results of magnetoconductivity measurements in Gain As quantum wells are presented. The 
observed magnetoconductivity appears due to the quantum interference, which lead to the weak 
localization effect. It is established that the details of the weak localization are controlled by the spin 
splitting of electron spectra. A theory is developed which takes into account both linear and cubic 
in electron wave vector terms in spin splitting, which arise due to the lack of inversion center in the 
crystal, as well as the linear terms which appear when the well itself is asymmetric. It is established 
that, unlike spin relaxation rate, contributions of different terms into magnetoconductivity are not 
additive. It is demonstrated that in the interval of electron densities under investigation ((0.98 — 
1.85) • lO^'^ cm~^) all three contribution are comparable and have to be taken into account to achieve 
a good agreement between the theory and experiment. The results obtained from comparison 
of the experiment and the theory have allowed us to determine what mechanisms dominate the 
spin relaxation in quantum wells and to improve the accuracy of determination of spin splitting 
parameters in A3B5 crystals and 2D structures. 
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I. INTRODUCTION 

The effect of the weak localization in metals and semi- 
conductors is caused by the interference of two electron 
waves which are scattered by the same centers (defects 
or impurities) but propagate in opposite directions along 
the same closed trajectory, and, therefore, return to the 
origin with equal phases. This effect increases the ef- 
fective scattering crossection,-, and, therefore, leads to 
a suppression of conductivitytl^d. In a magnetic field, 
the two waves propagating in the opposite directions ac- 
quire a phase difference 2e$/c, where $ is the magnetic 
flux through the area enclosed by the electron trajectory. 
This phase difference breaks the constructive interference 



and restores the conductivity to the value it would have 
without the quantum interference corrections. This is ob- 
served as an increase in conductivity with magnetic field, 
the effect known as positive rmuinetoconductivity (PMG) 
or negative magnetoresistancei'U. 

When spin effects are taken into account, the interfer- 
ence depends significantly on the total spin of the two 
electron waves. The singlet state with the total spin 
J = gives negative contribution to the conductivity 
(antilocalization effect). The triplet state with J = 1 
gives positive contribution to the conductivity. In the 
absence of spin relaxation the contribution of the singlet 
state is canceled by one of triplet states. As a result, 
the magnetic field dependence of the conductivity is the 
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same as for spinless particles. However, strong spin re- 
laxation can suppress the triplet state contribution with- 
out changing that of the singlet state, hence the total 
quantum correction may become positive. The interplay 
between negative magnetoconductivity at low fields and 
positive magnetoconductivity at high fields can lead to 
appearance of a minimum on the conductivity - magnetic 
field curve (antilocalization minimum). 

It was shown in the early papers by Hikami, Larkin, 
and NagaokaH and by Altshuler, Aronov, Larkin, and 
KhmelnitskiiJ that the behavior of the conductivity in 
weak magnetic fields depends essentially on the mecha- 
nism of the spin relaxation. Three mechanisms were con- 
sidered: Elliott- Yafet, otherwise known as skew scatter- 
ing mechanism, scattering on paramagnetic impurities, 
and Dyakonov-Perel mechanism, which arises from the 
spin splitting of carrier spectra in non-centrosymmetric 
media. 

Dyakonov-Perel mechanism is dominant in most A3B5 
cubic semiconductorsQ, with the exception of those with 
narrow band gap Eg and large spin-orbit splitting of 
the valence band A (for example, InSb). The same can 
be said about the low-dimensional structures fabricated 
from these materials. Presence of a antilocalization min- 
imum on the (j{B) curves in quantum 2D systems is a 
definite sign that the dominant spin-relaxatioii mecha- 
nism is the Dyakonov-Percl one. It is knowrfl that for 
the Elliott- Yafet mechanism in 2D structures the contri- 
bution of the singlet state with J = is exactly canceled 
by one of the triplet states, the one with J = 1 and 
J. = 0. 

Unlike bulk crystals, where the spin splitting is pro- 
portional to the cube of the wave vector A;, in 2D struc- 
tures the splitting has also terms linear in k (this is also 
true for strained A3B5 crystals and for hexagonal A2B6 
compounds). Furthermore, there are two linear in k con- 
tributions of essentially different nature. The first one, 
which arises from the lack of inversion in the original 
crystal (like the cubic term), is known as Dresselhaus 
terma, while the second, Rashba term, is caused by the 
asymmetry of the quantum well or heterojunction itsellEI. 

The direct measurements of the spin splitting using 
the Raman scattering in GaAs/AlGaAs quantum wellsEj 



have shown that, for electron densities Ns 
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both linear contributions are comparable. All three terms 
give additive contributions to the spin relaxation rate. 
When only the cubic in k term is present, its effect on the 
PMC is determined only by the spin relaxation rate, sim- 
ilarly to the two other spin relaxation mechanisms, and 
is described by the theory of Refs. In the presence 
of the linear in k term in the spin Hamiltonian it is nec- 
essary to take into account the correlations between tiie 
motion of electrons in the coordinate and spin spacestil. 
In the theory of coherent phenomena these correlations 
were first taken into account using the language of the 
spin-dependent vector potential in Refs. where 
this concept was appliadnte consideration of spin-orbit 
conductance oscillationgi3't3 and of the spin orbit effects 



in the universaLconductance fluctuation and persistent 
current in ringaiJ. In the theory of the anomalous magne- 
toconductivity the correlation between motion in real and 
spin spaces was first taken into account in Ref. It was 
shown that when linear and cubic in k terms are present, 
their contributions to spin phase-breaking are not addi- 
tive. Furthermore, as it was demonstrated in Ref. 16, the 



contributions of Rashba and Dresselhaus terms are also 
nonadditive, and the magnetoconductivity is determined 
not by their sum, but rather their difference. Similar 
effect, when spin-orbit phase-breaking may become neg- 
ligible due to the correlation of motion in real and spin 
spaces occurs in quasi- ID case and leads to two types 
of pronounced oscillations in the universal conductance 
fluctuations in ringsliZI. 

First experimental studies of the PMC in quantum 
2D structures were done in Refs. [l8|-^. Recently, the 
experimental observations of very pronounced effects of 
spin scattering-|(Mi weak localization conductivity correc- 
tions in GaA£2l'L3 and InAsE^ heterostructures have been 
reported. Different spin relaxation mechanisms are in- 
voked to explain the experimental data. In the Ref. ^ 
the spin relaxation is interpreted in the framework of 
the Dyakonov-Perel mechanism based on the bulk GaAs 
Hamiltonian. According to the authors, agreement with 
experimental data is achieved if one neglects in the spin- 
orbit Hamiltonian the term linear in the in-plane wave 
vector. This is in contradiction with the theoretical 
predictionst£l which show that, at least for low carrier 
concentrations, the linear term should be the dominant 
one. In Ref. ^ the same Dyakonov-Perel mechanism is 
used, but it is based on the Rashba term. In Ref. ^ it is 
assumed that the dominant mechanism of spin relaxation 
is the Elliott- Yafet scattering one. Recently, we have 
reported the measuremente-|Of the magnetoconductivity 
of GalnAs quantum wellsEJ. We have used Altshia.ler- 
Aronov-Larkin-Khmelnitskii (AALKh) calculationcl of 
quantum corrections to conductivity to interpret the ex- 
perimental data. However, the spin relaxation time, 
which enters the AALKh expressions, was calculated tak- 
ing into account not only linear but also cubic Dressel- 
haus terms of the spin splitting. We have demonstrated 
that, using this simplified theoretical approach, one can 
obtain right order of magnitude for the experimentally 
observed spin relaxation rates. 

In this work we present detailed experimental study 
of the negative magnetoconductivity in selectively doped 
GalnAs quantum wells with different 2D carrier densi- 
ties. We interpret them in the framework of recently de- 
veloped theory of the anomalous magnetoconductivi ty , va 
quantum wells which corrects the AALKh approachHE^. 
Comparison of experiment and theory allows to deter- 
mine importance of both Dresselhaus and Rashba terms 
for 2D systems. 
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II. THEORY 



(6) 



A. Spin relcixation 

The theory of the positive magnetoconductivity 
(PMC) for structures with spin sphtting hncar in wave 
vector was described very briefly in Refs. Below 
we present an outhne of this theory in more details. The 
spin splitting of the conduction band in cubiC[-|Crystals 
A3B5 is described by the following HamiltonianQ: 



i + 'i 



(1) 



where Ui are the Pauli matrices. In [001] quantum wells 
the size quantization gives rise to the terms in the Hamil- 
tonian Tis , which are linear in the in-plane wave, vector 
k ~ {kx,ky), in addition to the cubic termsHEa. The 
corresponding HamiltoniajB-for the conduction band elec- 
trons can be written aalaM 



(2) 



where cr — (ax,(yy), ^ — (^xi^y) are two-dimensional 
vectors with components in the plane of the quantum 
well. Vector 2Cl/h has the physical meaning of the pre- 
cession vector: its length equals the frequency of the spin 
precession and its direction defines the axis of the pre- 
cession. The spin splitting energy is equal to 2f2. To 
treat the spin relaxation problem in the case when is 
anisotropic in the 2D plane one has to decompose it into 
orthogonal spherical harmonics: 



fllx — COS (fi, fl: 



-f^a cos3(y5, 



fliy ~ ill siiKf, ilsj^ = — r23sin3(p, 



(3) 



f)« = Ik ^{kl) - Ifc^) , 



where fc^ = + ky, tanip — k^/ky, and (fc^) is the 
average squared wave vector in the direction z, normal to 
the quantum well (in this paper we take h = I everywhere 
except in final formulas). 

The spin salitting given by Eq. (^) represents the Dres- 
selhaus termO. In asymmetric quantum wells the Hamil- 
tonian Ti. contairualso terms of different symmetry, i.e. 
the Rashba termaa: 



The expressions for aQ and 7 are given in the Ap- 
pendix. The barriers of the well give rise to another 
contribution, usually also linear in £, which depends 
strongly on tha-details of the boundary conditions at the 
heterointerfacecSEZl. 

Both terms Eqs. (|) and (|) give additive contributions 
to the spin relaxation rate 1 / Ty , which is defined as 



'dt 



(7) 



where Si is an average projection of spin on the direction 
i. These contributions are 
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— = ^ = 2 {nln + nlrs) , 



(8) 



9 (1) (2) 

where flf — fl\ ' + fl\ and r„, n = 1, 3, is the relax- 
ation time of the respective component of the distribu- 
tion function fn{k) ^ cosn{ipf, +ipn) (V'n is an arbitrary 
phase) : 



— = / W((p) (1 — cosnip) dip. 



(9) 



Here W{^) is the probability of scattering by an angle 
If it does not depend on -d, all scattering times are equal 
to the elastic lifetime 



1 



To 



W{ip) dip. 



(10) 



When small- angle scattering dominates, 1 — cosmos 
{nipY /2 and 



{n > 1). 



(11) 



Formula Eq. (g) shows that the different harmonics of 
the precession vector add up in the spin relaxation rate 
with the weight equal to the relaxation times t„. Un- 
like the spin relaxation, the contributions of the different 
terms in the spin splitting into PMC are not additive. 



Furthermore. 



at n[^'^ 



±n[^^ and f]. 



the contri- 



butions of the two linear terms Eqs. (^ and (^ exactly 
cancel each other, and the magnetoconductivity looks as 
if there were no spin-orbit interaction at all. Analogous 
effect occur in weak localization conductance in wirescj. 



This term can be included in the Hamiltonian Eq. 
one includes additional terms into fl: 



(4) 
t) if 



— rj^^-* sin ip, n 



- _o(2) 



n\^> cos ip, ^ak. (5) 



In an uniform electric field £, the constant a is propor- 
tional to the field: 



B. Weak localization in two-dimensional structures 

The weak localization coatribution to the conductivity 
is given by the expressionBa: 

Aa^-^-2n,,orSj2 I Cc/3/3a(q) t^-^, (12) 
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where a and /3 are spin indices, g^^^, 
v'^Ti/2 is the diffusion coefficient, and = m/2iT is the 
density of states at the Fermi level at a given spin projec- 
tion. The matrix Gapi3a{(l) is called Cooperon and can 
be found from the following integral equation: 



Cai3is{k,k' ,q) = 



k.k' 



{Dti)~^, D = Therefore, the solution of inhomogeneous equation 
Eq. ( p!?!) will have large harmonic Cq, while the others 
will be small, because they appear due to presence of 
small terms in q and Q. Since the right-hand side of 
Eq. (^) contains linear and cubic in g terms, in is neces- 
sary to take into account only first and third harmonics. 
(13) From Eqs. @ - @ it follows that 



X Cxx'js{g,k',q). 

Here ^' is a scattering matrix element (including the 
concentration of scatterers) , which we assume here to be 
diagonal in spin indices. It is connected with W{(p) in 
Eq. (no) by the expression 



1^0 



Vi 



k.k' 



G^{Lj,k) are the Green's functions 



G^iu;, k) = {lu- E{k) - {aO.) ± — 

111 fc^ 

- = - + -, E{k)^—., 
Tf To 2m 



(14) 

(15) 
(16) 



Tip is the inelastic scattering time. After the integration 
by E{g) in the right-hand side of Eq. (|l^ the result is 
expanded up to second order terms in series in small pa- 
rameters tq/t^, vqtq, and Otq, where v = dE/dk. In the 
end, the following equation for the Cooperon is obtained: 



Vi 



k,k' 



27r 



Vi 



1 - i{vgq)To - i{(T + p)J7to - [Vgqfr'^ 



(17) 



-2(o-n)((pn)r2 - 2{vgq){a + p)nr^ ~ ^jc^ ^,(q). 

Here the Pauli matrices cr act on the first pair of spin 
indices a. A, while the matrices p- on the second pair 
/5A'. 

The equation 

Cfc^fc' = Aro 1 W{k',g)Cgf^,d^g (18) 
has the following harmonics as its eigcnfunctions: 

^k k' = '^(^fc ~ 'fi'^' ~ (19) 



According to Eq. ( |10[ ) the eigenfunction Co has the eigen- 
value Ao = 1, while other harmonics have eigenvalues 



A. 



1 



To 



(20) 



^g!k' = -*(^3 - ro){cT + p)J^3(g)C^ fc, 



(21) 



Here it is taken into account that there is a relation 
similar to Eq. ( [T^ ) for harmonics flia, (vq) ^ cos{ipg — 
ifq), and Qsa- 

Then we substitute C„ ^> = C° . , C^^^. , + , 
9,1^ g,K g^k g,k 

into Eq. (|^, and, using Eq. ( |2l| ) and retaining only the 

terms with zero harmonic, we obtain the equation for 

Co(9): 



1 



2'KVoT^ 



(22) 



where 



n = — 



1 1 



2 2 



{Vl\ti + ilaTa) (2 + o-^j/jj, -I- o-,yPy) 



(23) 



VTl 



f Py) 



In a magnetic field q become operators with the com- 
mutator 



where q± — qx^lsi iqy and 



5_ 



ieBD 

he 



(24) 



(25) 



This allows us to introduce creation and annihilation op- 
erators flt and a, respectively, for which [aa^] = 1: 

D'^\+^S'/\, D'^\.^S'/'a\ Dq'^6{aa^. 

(26) 

In the basis of the eigenfunction of the operator {aa^} = 
^(aa^ -fa^a) these operators have following non-zero ma- 
trix elements 



(n- l\a\n) = (?i|at |7i - 1) = 
{n\ {aat} \n) = « + ^• 



(27) 
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In a magnetic field, the integration over q should be 
replaced by summation over n. Then, 



Act 



(28) 



eigenf unctions of Jz (for rj^^-* = the basis is = 
(/i(n)F„+i,/o(n)F„,/_i(n)i^„_i)). Using the formula 



I A. 



D 



(36) 



where 

5 = 2^l^OTo' 51 C„/3;3a(n). (29) 

a,/3,n 

Since Eq. ( p2|) is essentially the Green function equa- 
tion its solution can be written as: 



C(n)"^ 



1 1 



where '^r,n and Er^n are the eigenfunctions and eigenval- 
ues of H: 



^r,n ^ r.n ■ 



(31) 



We now choose the basis consisting of the function 
^'o(a, /3), which is antisymmetric in spin indices and cor- 
responds to the total momentum J = 0, and of sym- 
metric functions which correspond to J = 1 and 
= m = — 1, 0, 1. According to Eq. (|30|), in this basis 
the sum in Eq. (Bq) is 



S 



where rim, 
ator Ti is 



and, therefore, 




1 



E 



1 



(32) 



Eo{n) ^ Ern{n) I 

^ ' rii— — l ' ' / 

1 /(5ri . For the term with J = the oper 
1 



Won = &{aa}} -I- 



(33) 



Ea{n) = 5 ( n + ij + -5-. 



(34) 



For the term with J = 1 we can use the relation Ji 
(cTi + pi)/2 to obtain 



n = 6{aa^} — + 2 {Vl\ti -f VlIt-^) (2 - J^) 



4zl]i^^r!i"Vi (4 - J2 



2((5ti 



,1/2 



(35) 



(2) 



where J± = (J^; ± iJy)/^/2. 

When f^f' = (or = 0), the operator 

Eq. (^ ) can be reduced to a block-diagonal form 
with 3x3 blocks if one uses the basis of functions 
= {fi{n)Fn-i, fo{n)Fn, f-i{n)Fn+i), where F„ are 
the eigenfunctions of the operator {aa^} and are the 



where \D\ is the determinant of 7i (Eq. (pSi)) and |-Dmm| 
are its minors of diagonal elements Dm^m the sum in 
Eq. (32) can be immediately calculatedEj. According to 
Eqs.ll), m, and 



AaiB) 



El 



/ 1 2ao + 1 + %^ 



B J " B 



(37) 



3a2 +2a„%i-l-2(2n+l 



so 
B 



{an + ^) a„-ia„+i - 2%^ [(2n + l)a„ - 1] , 



where C is the Euler's constant, 
1 iJ<^ i?50 

On = n -I- - + ■ ' 



2 B B ' 
B 



— ^Tfn, Hfr — 1 



^ '-"2, , on2_^ TJI _ t7-(i) u(2) 



(38) 



(1)^ 



r(2) 



2f^i'^ Tl, 



4feD 



and 5* is a digamma-function. 

If both fi^^^ = rip-* = and only the cubic in k term 
with 03 is present, the expression Eq. ( |37| ) can be further 
reduced to the formula, which was obtained earlier in 
Ref . 0: 



Acr(B) - Act(O) = 




1 -ff^ -H^so \ 

2 B B I 



B 



2 s 



(39) 



Note that, according to Ref. |T^, the value of Hso is twice 
that used in Ref. p. 

(1) (2) 

The case when Vl\ — ±£7]^ and f^a = is a special 



one. In this case the operator Eq. (23) is diagonal in the 
basis of functions ^'m if one uses coordinates x' \\ (110) 
and y' \\ (110): 



1 



D 



2 



{qy 



(40) 
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where <?y/„ = 2Q,i^JtiI Dm. Since the commutation re- 



lations Eq. (24) do not change when qyi is shifted by g^, 



the spin sphtting does not manifest itself in the magne- 
toconductivity, which is given by the simple formulaO 



Acr(B) - A(t(0) 



1 



(41) 



It was demonstrated in Ref. that this result appears 
because, when fl'"^^ = zbilj^'' and fl^ — 0, the total spin 
rotation for the motion along any closed trajectory is 
exactly zero. 

When Sl^^-* and rjj^-' are not equal or Q3 ^ 0, the only 
way to find eigenvalues is to diagonalize numerically 
the matrix 7i. The number of elements one has to take 
for a given value of magnetic field i3, or 5, is at least 
"^max = ^l^i'i and increases infinitely as B approaches 0. 
Note that the size of the matrix H is N = Sumax- For 
the detail of the numerical procedure, see Ref. HO. 



C. Elliott- Yafet spin relaxation mechanism 

It follows from Ref. ^ that in order to take into account 
the Elliott- Yafet spin relaxation mechanism one has to 
add a new term to the Hamiltonian Eq. ([3^): 



TisY = - 

where, according to Ref. 

1 1 



1 



T2 



2 J 



ZmEl (E, 



A 

g 3 



(42) 

(43) 
(44) 



and T2 is defined by Eq. (||). 

As a result, in the first and fourth terms of the for- 
mula for magnetoconductivity Eq. Hso should be 
replaced by Hso + Hey, where Hey is 
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FIG. 1. Sample structure (a) and band diagram (b) 
for GalnAs quantum well (sample Bl) as obtained from 
self-consistent calculations. The first two energy levels in the 
well are shown by solid lines, Fermi energy is shown by a 
dotted line. 



A. Samples 

Three AlGaAs/InGaAs/GaAs pseudomorphic quan- 
tum wells were studied. They were grown by the molec- 
ular beam epitaxy technique. The layer sequence of the 
structure was of the standard HEMT type and is shown 
in Fig. 1^. The two-dimensional electron gas was formed 
in the 13 nm thick InGaAs layer. Samples were 5-doped 
with Si (doping density Nd — 2.5 • 10^^ cm~^). Samples 
of the type A had a spacer thickness of 6 nm, samples of 
the type B had a 4 nm spacer and samples of the type C 
had a 2nm spacer. The samples had the Hall bar geom- 
etry with the length of 1.0 mm and the width of 0.1 mm 



H 



EY 



AheDTs 



It follows from Eq. (E^) that 



tltr T2 



III. EXPERIMENTAL PROCEDURES 
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TABLE I. Sample parameters: electron density Ns, mobil- 
ity /X, transport magnetic field Htr (Eq. (p^, and momentum 
relaxation time t\. 
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FIG. 2. Energies determining the Dresselhaus spin split- 
ting as a function of electron density Ns- Fermi energy Ef 
and quadrupled mean kinetic energy 4ii'z of the motion in the 
growth direction are shown by solid lines. Dotted line shows 
the difference 4_Bz - Ef that enters Eq. (||) for H^^^. 

with two current and four voltage probes. The distance 
between voltage probes was 0.3 mm. The samples were 
independently characterized by luminescence, high field 
transport, and cyclotron emission experimentaSa. The 
parameters are listed in Table |. In order to study the 
behavior of the structures as a function of electron den- 
sity TVs, the metastable properties of the DX-Si centers 
present in AlGaAs layer were employed. Different con- 
centrations were obtained by cooHng sample slowly in 
dark and then by illuminating it gradually by a light- 
emitting diode. This allowed us to tune carrier density 
from 0.98 • 10^^ cm~-^ to 1.95 • 10^^ cm^-^. We have mea- 
sured the Hall effect and Shubnikov-de-Haas oscillations 
to determine Ng and to verify that in all samples only the 
lowest subband is occupied. To calculate the energy lev- 
els in investigated quantum wells we first self-consistently 
calculate the 2_D wavefunctions, using the envelnpe func- 
tion approach in the Hartree approximationtJ'Ej. The 
potential entering into the zero-magnetic field Hamilto- 
nian takes into account the conduction band offset at 
each interface, and includes, in a self-consistent way, the 
electrostatic potential curvature due to the finite extent 
of the electron wavefunction. The boundary condition for 
the integration of Poisson equation within the 2D channel 
is the value of the built-in electric field in the buffer layer 
on the substrate side of the 2D channel. It originates 
from the pinning of the Fermi level near midgap in semi- 
insulating GaAs substrate. Any nonparabolicity effects 
on the effective masses were neglected. The calculations 
were performed for the temperature 4.2 K. Results of cal- 
culations are shown in Fig. With increasing concentra- 
tion both Fermi energy and kinetic energy of the motion 



in the growth direction increase. Their exact concentra- 
tion dependencies should be determined to calculate spin 
splitting and spin relaxation times. For every carrier den- 
sity N s the expectation value of the z-component of the 
kinetic energy was calculated. Figure ^ shows the result 
of such calculations for the quantum wells used in our ex- 
periments. We also show the Fermi energy as a function 
of carrier density Ng. 

B. Magnetic field generation and stability 

We have used a system of two superconducting coils 
(8T/8T) placed in the same cryostat. This system was 
earlier used to study cyclotron emission from the same 
sampl^A, B, and C and to determine their effective 
massesEil. The sample was placed in the center of the 
first coil. To generate the stable weak magnetic field, 
necessary for the antilocalization measurements we used 
a spread field of the second coil to compensate the field in 
the first one. The magnetic field scale was determined on 
the basis of measurements of the Hall voltages induced 
on the sample by both coils. Typically the constant mag- 
netic field in the sample coil was of the order of 400 Gauss 
and it was compensated by tuning the second coil field 
in the range from 12 to 14 kGauss. This way, both coils 
were operated in a stable and reproducible manner giv- 
ing in the sample space magnetic fields from - 30 Gauss 
to -1-30 Gauss. Small sample dimensions and the geome- 
try of the coils gave good magnetic field uniformity. We 
estimate thet the magnetic field have varied by less than 
0.1 Gauss over the sample. 

C. Conductivity measurements and temperature 
control 

We have used the standard direct current (DC) method 
to measure the conductivity with currents less than 20 
microampers to avoid sample heating. A high preci- 
sion voltmeter capable of measuring nV changes on mV 
signals was used to measure the conductivity and Hall 
voltages. The whole system was computer controlled. 
To avoid mechanical and temperature instabilities, the 
sample was not directly immersed in the liquid helium 
but was enclosed in the vacuum tight sample holder 
and cooled by helium exchange gas under the 50 mbar 
pressure. A calibrated Allan-Bradley resistor placed 
near the sample was used to measure the temperature 
which was stabilized between 4.2K and 4.3K. The exper- 
imental arrangement allowed simultaneous complemen- 
tary Shubnikov-de Haas and Hall effect measurements to 
determine carrier mobility and concentration for different 
sample illumination intensities. 
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FIG. 3. Experimental results (circles) and theoretical 
fits for the magnetoconductivity <j{B) — (j{0) for three dif- 
ferent samples: a) - Al, b) - Bl, and c) - C4. Solid lines 
show results of the theory outlined in Sec. |l[ Best fits ob- 
tained from Eqs. and are shown by dashed and dot- 
ted-dashed lines, respectively. Dotted vertical lines show the 
values B = 0.5Htr, which limit the intervals of applicability 
of all three theories. The fitting parameters H'^q, H'^q, Hso, 
and are given in Table 



IV. RESULTS AND DISCUSSION 



was developed in Rcfs. |3^-|35| for the case of isotropic scat- 
tering and with spin relaxation considered in the frame- 
work of AALKh theory. The comparison with the diffu- 
sion theory shows that in-magnetic field B = OAHtr the 
latter has an error of 6%ta. For the purpose of compar- 
ison with theory, we have selected only samples with B 
at the minimum of a smaller than OAHtr- 



A. General comments 

For all samples and for all carrier densities, the magne- 
toconductivity was a non-monotonic function of the mag- 
netic field. As we have mentioned before, presence of a 
minimum on the cr(B) curves is a definite sign that the 
dominant spin-relaxation mechanism is the Dyakonov- 
Perel one. For the Elliott- Yafet mechanism in 2D struc- 
tures, the contribution of the singlet state with J = is 
exactly canceled by one of the triplet states, namely the 
one with J = 1 and = 0, which is immediately evident 
from Eq. (^). Using Eq. ( ^6|) one can show that even 
for the highest density Ng — 2 - 10^^ cm~^ and ti/t2 = 4, 
the characteristic magnetic field Hey does not exceed 
4 • 10~^ Htr, which is much smaller than Hso- For the 
scattering on paramagnetic impurities, the negative mag- 
netoconductivity at lowest fields does not exist both in 
2D and 31? systems. 

As we have already noted, the theory presented in this 
paper uses the diffusion approximation, which is valid 
only when all of the fields H^ and Hso are smaller than 
Htr- Kinetic theory, which is free from this limitation, 



B. Description of fitting procedure 

The experimental data for each sample are fitted with 
the results of three differept theoretical models. First 
one is the AALKh theoryQ Eq. ( ^9| ) and has Hso and 
H^ as fitting parameters. The second one corresponds to 
the physical situation where one of the linear terms H'qq 
or H^Q dominates, and the Eq. (|37| ) can 1^ used with 



the fitting parameters Hso, H'^q, and Hji3- The last 
theory takes into account all the terms Hsoi h'^so^ ^^'^ 

(2) 

Hso exactly. The results of this theory were obtained 
by numerical di agona lization of the matrix Eq. (|3^), as 
described in Sec. [I B and Ref. nq. The fitting parameters 
in this case are Hso, H^so, nfl, and H^ (see Table 0). 



The fitting of the experimental data by Eqs. (^7|) and 
( |39|) was done by weighted explicit orthogonal distance 
regression using the software package 0DRPACKE3. The 
weights were selected to increase the importance of the 
low-field part of the magnetoconductivity curve. The 
calculation of the magnetoconductivity by numerical di- 
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agon aliza tion of the matrix Eq. (pq), as described in 
Sec. [IB , requires large amounts of computer time, and 
we could not afford to use the automated fitting with 
these results. The fitting was done "by hand" , using 
empirically gained knowledge on how changing different 
fitting parameters affect the magnetoconductivity curve. 



C. Experimental results 

In Fig. H a-c we show the results of the measurements 
of the conductivity a as a function of magnetic field for 
three different samples. To compare the results for dif- 
ferent carrier densities we plot cr(B) — a(0) in units of 
e^/2TT'^h = 1.2310"^ r^-i. The value of cr(B)-CT(0) gives 
the conductivity change induced by the applied magnetic 
field and can be directly compared with theory. The cir- 
cles show the experimental data, the results of the theory 
presented in Sec. || are shown by solid lines. The values 

of parameters Hso, ^so' ^"^"^ ^ ^^^^ values of 

Ns and Htr, are given in Table ||. 

Before the quantitative analysis of the experimental 
data, we would like to point out some of their general 
features. The position of the characteristic conductiv- 
ity minimum which shifts from 2.5 Gs in Fig. |^a to 5 
Gs in Fig. ^ is largely determined by the value of Hso, 
and, hence, by the spin relaxation rate. With increasing 
carrier density this minimum shifts towards higher 
magnetic fields. This indicates an increase in the effi- 
ciency of the spin relaxation. One can also observe that 
the minimum becomes more pronounced when the ratio 
Hso /Hip increases: the minimal value of a{B) — a{0) 
is about for the sample Al, 0.0le^/2n^h for 

the sample C4, but increases to O.IIe^/27r^?i for Bl. This 
shows that the magnitude of the antilocahzation effect 
depends strongly on the ratio of the phase-breaking and 
spin relaxation rates. Small phase-breaking rate and fast 
spin relaxation increase the magnitude of the antilocal- 
ization phenomenon. When the two rates are compara- 
ble, the antilocahzation minimum almost vanishes (this 
can be seen in Fig. ||c for the sample C4). 

In Fig. I^a for the sample Al the dashed line shows the 

~ i.e. with H^sl = 0. The 
-^so ~ O.OSGs is also close to 



best fit obtained using Eq. (^ 
best fit value of H'g^ - "^'^^ 
0. Hence, the dashed curve almost coincides with the 
dashed-dotted line, which shows the result of AALKh 
theory, Eq. ( |39| ) . Both theories fit the experimental data 
seemingly quite well. However, the values of parame- 
ters required to achieve this agreement {Hso ~ 0.8 and 

(2) 

Hso ~ 0) '^^ ^ sharp contradiction with theoretical 
calculations of Hso and experimental measurements of 
7, while the theory presented in this paper fits the ex- 
periment using the parameters a and 7 which agree with 
other measurements and calculations (see Sec. IV D , the 
Appendix, and Refs. |,|7|-|o|). 

The results for the sample Bl are shown in Fig. ^ b. 
Again, the dashed line shows the fit by Eq. ( p7| ) and the 



dashed-dotted line - by Eq. ( p9| ) . One can see that in this 

case the theory with both H^q and -ff^Q, presented in 
this paper (solid line), gives somewhat better agreement 
with the experiment in the vicinity of the conductivity 
minimum. The general agreement of all curves with ex- 
periment is of similar quality, but again in order to bring 
Eqs. (^7|) and ( ^Of ) in agreement with experiment one has 
to use unrealistic values of Hso and H'^q . 

Fig. II c shows the results for the sample C4. The 
dotted-dashed line in Fig. || c shows the result of AALKh 
theory, Eq. (|39|). One can see that for B > 10 Gs this 
curve deviates from the experimental results quite signif- 
icantly. For this sample, as well as for two other samples 



C2 and C3 with large electron densities and H 
we have taken to be equal to its theoretical value for 
= 24 eV A"^. One can see from Fig. ^ c that the solid 



7 



curve, computed for -ff^Q 
practically coincides with the curve, computed using 
Eq. (H) for H^gl = and H^g^ = 3.97 Gs. This means 
that for large Ns the experiment allows to measure only 



0.34 Gs and H^s^l = 4.32 Gs, 



the difference H^q — H^so- discussion above shows 
that the new theoretical approaches developed in this 
work allow to improve the description of the magnetocon- 
ductivity dependencies and to obtain meaningful param- 
eters from the fits. In the next section we show that us- 
ing the compkete theoretical description with H^so^ ^^SO^ 
and Hso as the parameters one can get a consistent de- 
scription of experimental data for samples with different 
carrier densities. 



r(i) 



D. Carrier density dependencies 

In Fig. ^ we show the values of H^so ^'^'^ ^SO ^ 
function of iV^ for all samples we have studied, as ob- 
tained from the fitting of the experimental results by our 
theory. We also show the theoretical curves for these 
fields, calculated using Eqs. (||-|) and (H): 



H 



(1) 
so 



H 



(2) 

SO 



V2alN^ 



m 


r 




- 1 


mo 














+ 1 



(47) 



where A'o is the charge density in the depletion layer, 
Ez = Ti'^{k1)/2m is the kinetic energy of motion in z- 
direction, and 



n crriQ 
Aeh^ 



(48) 



Here too is a free electron mass. The calculations are 
done for 7 = 24 eVA'^ and ao = 7.2 A^. These values 
allow a good description of the experimental data and 
are close to those obtained from k ■ p and tight-binding 
calculations for Gao.85Ino.15As (see the Appendix). The 
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TABLE II. Parameters of the best fits for three samples Al, Bl, and C4 (shown in Fig. |3| a), b), and c), respectively) as 
obtained from the theory of Sec. ^ (rows I), from the theory of Ref. |l^ and Eq. ^ (rows II), and from the theory of Ref. q and 
Eq. (rows III). All magnetic fields are in Gauss. 



Sample 


Theory 






Hso 




H's^ = Hso ~ Hl.>, - Hl.>, 




I 


0.62 


1.41 


2.69 


0.66 


0.66 


Al 


II 





0.03 


0.85 


0.66 


0.82 




III 








0.77 


0.59 


0.77 




I 


0.66 


1.91 


3.52 


0.60 


0.96 


Bl 


II 





0.87 


1.89 


0.58 


1.02 




III 








1.08 


0.53 


1.08 




I 


0.34 


4.32 


5.98 


3.03 


1.33 


C4 


II 





3.97 


5.30 


3.03 


1.51 




III 








2.18 


2.38 


2.18 



ratio Ez/Ep is calculated using Fig. 0. When calculat- 

(2) 

ing HgQ, we have assumed that the average field in the 
well is one half of the maximum field £ = AireNs/K. We 
have also taken into account the charge in the depletion 
layer A^o = 0.58 • _10^^ cm~^. The value of ao was calcu- 



lated using Eq. (Al). If one takes into account the barri- 
ers, using theory of Refs. and the self-consistently 
calculated wave functions, the value of ao will increase 
by about 60% for the electron densities in the interval 
Ns — (1 — 2) ■ 10^^ cm~^. This would increase the value 

(2) 

of HgQ approximately 2.5 times, but such large values 

of HgQ clearly do not agree with the experiment. It is 
likely that the barrier contribution depends very strongly 
on their microscopic structure, which may be very differ- 
ent from the abrupt interface model, used in the theory. 
It is also plausible that the different barrier structure is 

(2) 

responsible for the relatively large value of HgQ for the 
sample CI, for which ao = 8.8 A. 

It can be shown using the Eq. (|^) and data of Fig. || 
that for Ns < Nso = 7 ■ 10^^ cm^^ the Dresselhaus term 
decreases with increasing Ng, vanishes for Ns ~ Nso, 
and then begins to increase. One can see from Fig. ^ 
that for A'^s > 1 • 10^^ cm~^ the Rashba term exceeds the 
Dresselhaus term. Consequently, we denote the larger 



contribution in Fig. |4| as ff *"^^ 



One can see from Figure W that the general character 

(1) (2) 

of the density dependence of Hg^ and H^q agrees with 
the theory, and their values are close to those calculated 
using the above values of 7 and ao. 

In Fig. ^ b we show a similar density dependence but 
for the cubic in k Dresselhaus term Hso — t^so- 
theoretical formula for this field is 



Hso — -ff = 



so 



mi 



m 
Too 



N: 



(49) 



The top curve corresponds to T1/T3 — 1 and the bottom 
one - to ti/ts = 2. 

In the case of isotropic scattering, which is the case of 
short range potentials scattering, probabillity W{ip) in 
formula Eq. (|^) is angle independent and ti/ts = 1. If 



only small angle scattering is important (that is the case 
of scattering by the Coulomb potential) then ti/ts = 9 
(see Eq. (pT|)). In our case we find ti/t^ to be in the 
range from 1 to 2. It is probably because scattering 
in our samples is the mixture of short and long range 
scattering. The short range scattering is probably due 
to alloy scattering that is known to be mobility limiting 
mechanism in GalnAs quantum wells. Long range scat- 
tering is most probably due to scattering on the ionized 
impurities in the (5-doped layer. Role of scattering by the 
charged impurities in the i5-doped layer was confirmed by 
observations of charge correlation effects (see Ref. |2l) . 



V. CONCLUSION 

In conclusion, we have presented new experimental 
studies of positive magnetoconductivity caused by the 
weak localization in selectively doped Gain As quan- 
tum wells with different carrier densities. The com- 
plete interpretation of the observations is obtained in the 
framework of recently developed comprehensive theory 
of quantum corrections to conductivity. In this theory, 
we correctly take into account both linear and cubic in 
the wave vector terms of the spin splitting Hamiltonian. 
These terms arise due to the lack of the inversion sym- 
metry of the crystal. We also include the linear splitting 
terms which appear when the quantum well itself is not 
symmetric. 

It is shown that in the density range where all the 
above terms are comparable, the new theory allows not 
only to achieve good agreement with the experiment but, 
unlike earlier theories, also gives the values for the param- 
eters of the spin splitting which are in agreement with 
previous optical experimentsErliy and theoretical calcula- 
tions. Therefore, our research answers the question what 
spin relaxation mechanism dominates for different elec- 
tron densities and how it should be taken into account to 
describe the weak localization and antylocalization phe- 
nomena in quantum wells 



10 



TABLE III. Values of the parameters for GaAs and InAs calculated using the sp^s* model and the results of the k ■ p 
model. The parameters oi k-p model were taken from Ref. except those marked by asterisk, which were taken from Ref. ^ 
Parameters for Gao.85Ino.15As were obtained by linear interpolation of the k • p model parameters between GaAs and InAs. 
The values of 7 and ao as obtained in these models are also given. 





n A 




Tn A 11. 






A 
I5AS 




K ■ p 


Sp-iS 


fe ■ p 


spss 


K • p 


SpsS* 


J^g (eV) 


l.Oiy 




n /1 
U.4z 


U.4if5 


1.00 


i.004 


A (eV) 


0.341 


0.341 


0.38 


0.38 


0.347 


0.347 


K (eV) 


2.97 


2.98 


3.97 


3.95 


3.12 


3.104 


A' (eV) 


0.171 


0.159 


0.24 


0.26 


0.181 


0.20 


P (eV A) 


10.49* 


10.23 


9.2* 


9.22 


10.29 


10.16 


P' (eV A) 


4.78* 


1.46 


0.87* 


1.06 


4.20 


1.03 


(eV A) 


-8.16* 


-7.0 


-8.33* 


-7.27 


-8.18 


-7.03 


m 
mo 

7 (eVA^) 


0.0665 


0.066 


0.023 


0.023 


0.06 


0.06 


27.5 


10 


26.9 


71 


27.7 


13 


ao (A^) 


5.33 


5.15 


116.74 


118.5 


7.2 


7.05 



^The sign of Q in fc ■ p model is not defined can be chosen to be sthe same as in spss* model. 
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APPENDIX A: SPIN SPLITTING IN GAAS, 
INAS, AND GAINAS 

Below we present results of the calculations of ao and 
7 for GaAs, InAs and Gao.85Ino.15As in k ■ p and tight- 
binding calculations. The tight-binding calculations were 
done in the 20 baiid|rtight-binding model including the 
spin-orbit couplingEZIEa. Our calculations of electronic 
properties use sp^s* tight-binding parameters especially 
chosen so as to reproduce several features of the funda- 
mental properties of bulk constituents. We state some 
analytical relations connecting the effective masses and 
the deformation potentials at the F point, and the fifteen 
parameters of the spss* 20 band tight-binding modeled. 
Using these relations, as well as other relations bntween 
the fifteen parameters and F and X energy valueg23, we 
get a set of parameters which accurately reproduces the 
effective masses at the F point, the [001] deformation 



potential and overall band structure in accordance with 
refiectivity and photoemission measurements^. 

Such a procedure has been already checked to give a 
good description of refiectivity data-in uniaxially stressed 
GaAs/Gao.sglno.ii As superlatticeal3. In this work we use 
it to obtain InAs and Gao.85Ino.15As parameters. Us- 
ing these parameters, we calculate ao on the basis of 
Eq. (Al). In order to determine the value of 7 we calcu- 
late the value of the spin splitting as a function of k along 
(110) direction. It follows a cubic dependence in k from 
which we extract the values of 7 given in Table FT^. T he 
parameters of the k-p model were taken from Refs7||,||. 
The k p parameters for Gao.85Ino.15As were obtained by 
linear interpolation between GaAs and InAs. 

In the 3-band k ■ p model one takes into account the 
states of the conduction band Fi ( Fg) with the Bloch 
function S, the valence band Fi5„ (Fs -I- F7) with func- 
tions X, Y, Z, and the higher band Fi5c (Tsc + ^Tc) with 
functions X', Y', and Z' . The energies of these states at 
k ^ are: Er, = 0, Er, = -Eg, Er, - -{Eg + A), 
E^^^ = Eg, and Spsc — + ^'- ^^^^ ^^'VlV^/^/Ptr^' 
and ao are given by the following expressionaH'CJ aOE3: 



mo 



2 mo] SEgJ- 2A ^ -.^g 



3E' + A' 



3 [' EgiEg + A) E'g{E'g + A') 

4 PP'Q / A A' 
^^'3Eg{E'g + A') [E;TA^Wg 

ao = l[p' [E;' - {Eg + A)-'] - 



(Al) 



P' 



where P = ih/mo {S\p^\Z) , P' = ih/nio {S\p,\Z') , and 
Q — ih/niQ {X\py\Z') are the interband matrix elements, 
mo is the free electron mass, p = —ihX7. Here we do not 
take into account the contribution into 7 and ao which 
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FIG. 4. Characteristic magnetic fields as a function of the 
electron density Na- 

a) - Density dependencies of the Dresselhaus (H^gQ) and 
Rashba {H^q) linear terms are shown by dotted and solid 
lines, respectively. Calculations were done according to Eq. ^ 
with 7 = 24 eVA^ and oo = 7.3 A^. The values of these fields 
as obtained from best fit with the Sec. ^theory are shown by 
squares (Rashba term) and circles (Dresselhaus term). 



(1) 



(2) 



b) - The Dresselhaus cubic term Hso — H^q — J/gT, as a 
function of Ns. The lines are calculated using Eq. pq for 
7 = 24 eVA^. Solid line shows results for an isotropic scat- 
tering, ti/t3 = 1, Dotted line - for ti/t3 = 2. 



arises from spin-orbit mixing of the states Ti^y and Ti^c- 
The values of 7 obtained for GaAs from tight-binding 
calculations are usually smaller then those given by the 
k ■ p model (see Table For example, in Ref. |3^ the 
tight-binding calculations give the value 27 = 17.8 eVA^. 
In the later work Ref. ^ for the same parameter 27 (this 
time called 7) authors obtain the value 17 eVA^. Cal- 
culated values of 7 should be compared to the experi- 
mental values 24 eVA'^ for bulk GaAsD and to the if^ 
cently obtained value for GaAs/GaAlAs quantum wellsEj 
16.5±3 eVA^. 
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